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TECHNICAL MEMORANDUM 78220 


PROPELLANT GRAIN DYNAMICS IN AFT 
ATTACH RING OF SHUTTIE SOLID 
ROCKET BOOSTIR 

I. SUMMARY 

Tho dynamic modulua of solid propellant grains Is known to depend 
significantly on temperature and frequency, but equally Important Is the dynamic 
strain amplitude In excess of lO”**. This can be a troublesome property In 
dynamics analysis because the absolute amplitude of motion throughout the 
grain structure must be anticipated In order to select the real modulus from 
propellant test data. 

A practical method Is presented for determining the Solid Rocket 
Booster (SRB) grain modulus In the aft attach ring of Shuttle-mated vehicles 
during lift-off. It Is essentially a plane unsymmetrlcal ring mode analysis 
using an existing structural dynamics program, but Inherent procedures may be 
extended to other grain vlbi’atlon problems. 

The analytical method Is demonstx'ated In the resolution of the flight 
vehicle dynamic modulus over a probable range of pitch excitations during 
lift-off. 


II. INTRODUCTION 

The unusual tangential mounting of the SRB onto the external tank ( ET) 
at the aft attach ring was suspected of causing such large deformations of the 
grain cross section during pitch excitation that a critical selection of the pro- 
pellant stiffness was required for the d 3 'namlc math model. 

The purpose of tills study was to specifically determine the dynamic 
real modulus existing in the SRB flight vehicles during lift-off. The analysis 
was limited to the cross section of the grain In the ring but inherent procedures 
would be extended to other grain structures If Indicated. 



The predominant effort was In developing practical techniques to deter- 
mine the dynamic strain of the grain as a function of absolute modal amplitudes 
which had to bo related to external excitation. Additional propellant tests were 
required to Include frequencies from 1 to 5 IIz with a range at least an order 
of magnitude of dynamic strain variations and at different temperatures. 

The Mated Vehicle Ground Vertical Test (MVGVT) data were used to 
determine the actual behavior of the SRB aft attach ring. The method for 
determining the dynamic strain and dynamic modulus was demonstrated using 
the excitation data at the ring and the fundamental frequency of the mated sys- 
tem. Ultimately, the MVGVT data and live propellant properties at extreme 
launch temperatures were projected to determine the grain modulus at the 
attach ring on flight vehicles during lift-off. 

It Is hoped that the liberal documentation of this study will provide 
useful Insights for application to future propellant Improvement programs. 

III. FREEVIBRATION OF RING-GRAIN COMPOSITE STRUCTURE 


The propellant grain outer diameter Is bonded to the Inner surface of 
the SRB aft attach ring. The ring Is tangentially supported and accelerated as 
sketched In Figure 1. A free vibration analysis of this composite structure 
was conducted to Identify and define relationships of all pertinent dynamic 
parameters affecting the grain stiffness before determining the dynamic modu- 
lus existing during forced excitation. 

This analysis Included a screening of environments affecting dynamic 
parameters, a development of absolute dynamic strain from modal analysis, 
and testing considerations for deriving supplementary propellant data. This 
part of the study progressed In that order and only a minimum of details are 
excluded. 


A. Dominant Operational Environments 

The operational environments that have the most effect on the propellant 
grain stiffness in the aft attach ring are those existing during lift-off; i. e. , 
grain temperature, pitch acceleration, and combustion pressure. 
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• EXTIfMAL ACCCLiRATION 



KlRuru 1. (train cross section displacement due to Z-a\is 
acceleration at SHB rin^ station 1511. 


The propellant real modulus increases slf^nlficantly over the anticipated 
ambient range of 90* to 40* F. For any given lift-off day, the grain temperature 
is reasonably uniform throughout the cross section so that only dynamic proper- 
ties of the propellant at constant temperatures need to be considered. 

The combustion pressure aiding in t.he grain cavity imposes no direct 
grain strain. Propellant has a Poisson ratio of nearly 0. 5 and is restrained 
by the relatively stiff metal ring so that no grain displacements can occur due 
to pressure. However, the incompressible grain will transmit the coml)Ustior. 
pressure faithfully to the laterally vibrating ring which reduces the amplitude 
through nonlinear shell effeids. This reduction of ring vibration amplitude 


reduces the grain strains to loss than tliat of no combustion pressure. This 
phenomenon Is accommodated by the modal analysis program. 

Tho third significant environment Is the pitch excitation during lift-off, 
which Is caused by the sudden release Of the vehicle from on on-pad sprung 
condition due to winds and main engine thiaist. The excitation component of 
interest Is tlie acceleration acting along tho Z-axls which, for tho purpose of 
free vibration analysis. Is conceived as an external force that defines tho 
Initial displacement of natural vibrating ring-grain structures. The magnitude 
of this Z-axls excitation Is related to the grain dynamic strain. 

Each of these thx’oe parameters is assumed constant for a configuration. 
A change In any one or all parameters establishes a new configuration. 


B. Representative Dynamic Strain 

The NASTRAN modal analysis of tho grain cross section bonded to the 
aft attach ring was modeled with finite elements Illustrated In Figure 2. For 
a given set of Input conditions, 1. o. , combustion pressure and propellant shear 
real modulus, G* , the program calculated the natural mode frequency and rela- 
tive amplitudes of each finite element node In the Y and Z axes noted as T2 and 
T3 eigenvectors, respectively. The program was also used to determine the 
absolute amplitudes of element nodes for the statically displaced ring-grain 
structure at one earth gravity acceleration. Both established the method for 
determining an absolute representative dynamic strain of tho propellant grain, 

1. Relative Strains from Modal Analysis. A typical finite element with 
notations and conversions used to define radial, tangential, and shear strains 
Is presented In Figure 3. Notations and symbols used In Equations (l) through 
(8), though not conventional to structural analysis, were devised to conform 
to Tektronix computer programming. Because of the repetitive nature of 
applying these equations to each finite element and averaging, this program was 
developed and a printout is provided in the appendix. 
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POSITIUM NliMSfR 0> T2& TI 
[lOINvrCTORS PROM MATH MODf L 
POR fCCTlQN bOTOO 



rijjuro 2. SHH urain cross suction with finite elemont 
identification and typical strains. 


TANGENTIAL (MSALACEMENT 
M - nsiN^> TlCOt^ 




Figure 3. Finite elemunt (typical) with coordinates and notations 
used in Equations (1) throuKh (8). 


fi 



Convoi'slon of olgonvoctors to rolattvo radial and tangential displacements, 
L and M respectively, was obtained from 


LX = T2(J + 1, 1+ 1) cos^(l) -H T3(J +1, I+l) sin 0(1) 
L2»T2(J+1, 1+2) cos 0(2)+ T3(J+X, 1+2) sin 0(2) 
L3 => T2(J, 1+ X) COD 0(X) + T3(J, I +X) sin 0(1) 

L4 = T2(J, 1+2) cos 0(2) + T3(J, 1+2) sin 0(2) 








(la) 


MX « T2(J+ X, I +X) Sin 0(X) - T3(J+ 1, 1+ l) cos 0(X) 
M2= T2(J+ X, 1+2) sin 0(2) - T3(j+ X, 2) cos 0(2) 
M3= T2(J, 1+ X) sin 0(1) - T3(J, 1+ l) cos 0(l) f 


. (lb) 


M4“ T2(J, 1+ 2) sin 0(2) - T3(j, 1+ 2) cos 0(2) 


The relative radial strain of each element was defined by 


UX~ 


( LX + L2 >■ L3 - L4) 
2S ( J) 


( 2 ) 


the relative tangential strain was defined by 


__ (Ml - M2) (M3-M4) 

2B(J+1) 2B(J) 


( 3 ) 


and the relative shear strain was defined by 


( Ml + M2 - M3 - M4) ( LX - L2) ( L3 - I/l) 

2S(J) ‘^2B(J+X)'’' 2B(J) 


(4) 
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whoro 


A(J) « [R2(J+ 1) 

Is tlio olomont area and 

S(J) « R(J+ 1) - R(J) 

B(J)» R(J) 27t/N 
B(J+ 1) « r(J+ 1) 27t/N 

are side lengths of the 'Jlomonts. The number of elements, N, dividing tlie 
cx’oss section circle Is IG. 

The three relative strains of Equations (2), (3), and (4), acting 
simultaneously, were combined Into a single familiar parameter (relative 
dynamic strain) tlirough the strain energy density concept from elastic theory. 
The strain energy per unit volume may bo expressed without stress terms by 

W= 0.6\e2+ + g 2^G+ 0.6y2 G» , 

whoro G’ is the shear modulus. The volumetric change, e, of a relatively 

incompressible propellant Is zero which eliminates the first toi*m of the equation. 

The shear strain, -y, can bo reduced to zero by replacing tlio radial strain, c , 

r 

and tangential strain, c^, with the maximum and minimum principal strains 

0.6 

FI, F3 = 0. 5 (m + U2) ± 0. 5 [(Ul + U2)2 + (U3)^] . (?) 

These are relative principal strains since relative strains of Equations (2), 

(3), and (4) are used. 
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Substituting tliosQ principal strains into tho strain energy density equation 
and rearranging, tlio relative dynamic strain squared for each element is given 

t)y 


(F4)2« (Fl)‘‘+ (F3)2=W/G' 

Tho relative dynamic strain of each finite element varies over tho 
unsymmotrically displaced grain cross sectloiit However, to reduce tho com- 
putational complexity, grain properties are averaged over tho grain total area. 
Accordingly, tho average relative dynamic strain Is 


M N 

2 2 [F1* + F22](I, J)A(I, J) 

J°1 I-l 

N 

N X A(J) 

J=4 


1/2 


( 8 ) 


resulting from modal analysis wltli a real shear modulus input of G' . 

2. Absolute Amplitude of Dynamic Sti’aln . If the absolute displacement 
of one finite element node point could bo established. It would tlion follow that 
the dynamic strain absolute value would equal Equation (8) by a proportionality 
factor. A method for obtaining the real displacement of a finite element node 
was realized by using tho same modal analysis program to obtain the absolute 
deflection of the supported ring-grain structure Illustrated In Figure 1 with a 
1-g static load acting along tho Z-axls, This Is a fair simulation of a known 
force Imposed on the structure to produce an absolute Initial displacement 
before releasing It to vibrato freely with that amplitude. It Is also a good 
representation of the ring force excitation at the attach fitting by a known pitch 
acceleration g during lift-off. 

The relationships of parameters between the static and modal system 
were derived from the normal mode equation of motion without damping 



I 

! 


whoro 7 ]j aro the normal eoordlnatae defined by the coordinate transformation 

w(y, t) = 2 0.(y)il.(t) 

l«l 


In which w(yf t) Is the deflection of a point In the displacement configuration 
cjb^(y) of a vibrating system In the l*-^ natural mode. Applying an external force 


]7(y t) = m(y , t) g and using deflections of the coordinate transformation 

J J 

equation, the virtual work Is 


m 

6W= 2 F(y , t) 5w(y,» t) 5 


and In normal coordinates it is 


m 

6W = g 2 ‘"(yJ 
j=l ^ 



0l(yj) 


Rearranging and equating with the virtual work obtained from the generalized 
force 


n n m 

6W= J q 6i) = g 2 6l, E '^’ 1 (^ 1 ) 

1=1 ‘ ‘ 1=1 ‘i=i ^ ‘ J 


one obtains 


m 

Q, = g E 

j=l ^ ^ 
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Substituting this Into tlie preceding equation of motion and letting = 0 
for static deflection and for generalized mass normalized to unity, the 

normal coordinate for tlio deflected system Is 


m 


-irf 

I j„l j j 


Identifying the static system by subscript *’0" and the fundamental mode 
system by no subscript, one can establish the ratio between the normal coor- 
dinates of the two systems as 


77 =T] 


0 



The static absolute displacement AZ^ of a finite element node point is 

related to the eigenvector T3 of the corrospondrng node of the fundamental 
mode system by 


AZo = 77^T3 


The node point used in this study Is identified by position No. 50209 in the 
modal analysis printout. The stiffness, k, is tlie combined metal ring stiff- 
ness, C, In parallel v/lth the propellant grain real modulus related by 

k=(C+G’) 

Substituting both expressions into the preceding normal coordinates equation, 
one obtains 



1 

0 

+ 

0 

0 - 
1 

r s 


. 

[c + G> 


(9) 
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which Is the desired factor by which to multiply Equation (8) to convert rela- 
tive strain to absolute dynamic strain 


e^ = T7(F4) . (10) 

Since all propellant properties are determined from shear tests, tlie 
tensile dynamic strain, Is converted to shear dynamic strain, by 

equating the tension strain energy density ofi](F4) In Equation (lO) with the 
shear strain energy density of a lap shear specimen 

W = t?^(F4)2 G» G*/2 


The resulting conversion Is 


=3 ^/T Tj ( F4) = ^/^ 

Solving forrj and substituting Into equation (9), the dynamic shear strain 
of the free vibrating grain structure with Initial condition g Is 


y^ = gDP 

Where the static parameters are lumped Into the constant 


( 11 ) 


AZ 


( 12 ) 


and g^ Is unity. The modal system parameters are expressed by 


P = 


F4 

T3(C + G») 


(13) 
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The metal ring stiffness In Equations (X2) and (X3) Is determined from 
two static deflection analyses. One uses only the ring with the propellant modulus 
equal to zero; the otlier uses some arbitrary value of Gjj greater than 300 to 

obtain AZ^ and AZ^ deflections, respectively, at the same finite element node 

point used In Equation (o). Equating the two results, the stiffness constant is 
calculated by 




(14) 


This constant must be calculated for each SRB ring configuration having a 
different stiffness due to changes In geometry or combustion pressure. 


C. Digest of Stiffness Parameters 

It Is now appropriate to briefly recall sources and models already dis- 
cussed by which the eight parameters Identified act simultaneously to define 
the propellant grain stiffness. The modal analysis must satisfy the environment 
conditions and be compatible with propellant characteristics. 

Three of tliese parameters originated from lift-off environments. The 
combustion pressure, p, Is an input to tlie modal analysis which effectively 
stiffens the metal ring. It Is also used in tlie static analysis to determine tlie 
ring parameter, C, by Equation (14). The pitch excitation, g, is assumed to 
initiate the free vibration absolute amplitude by which absolute dynamic 
strain, y, Is calculated In Equation (11). The ambient temperature, T, Is used 
to estimate the grain temperature which grossly affects the propellant stiffness 
as determined from propellant dynamic test. Results from any grain modulus 
analysis at lift-off must address these three environment parameters. 

The shear real modulus, G’ , Is tlie desired value to define the propellant 
grain stiffness In this study. It Is also the most significant propellant parameter 
Input to the modal analysis from which the fundamental frequency, oj, Is obtained 
as well as the eigenvectors by which the relative dynamic strain of Equation (8) 
Is calculated. Another real modulus, Gj^, is used In the static analysis to 

determine the ring constant, C, of Equation (14). 
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Of these eight parameters, tliroe are accounted for la the modal pro- 
gram and five are embodied In Equation (ll) but, temperature Is not directly 
required by either one. 

Finally, the restraining model of the four parameters already cited Is 
provided by the propellant behavior determined from dynamic tost data. These 
parameters are the real modulus, dynamic strain, frequency, and temperature 
(sometimes referred to as propellant properties) . They are related by 


G» = 


g (j 
„b e 


( 15 ) 


where a, b, c, and e are determined by the tost. However, these parameters 
are more conveniently plotted. 


IV. PROPELLANT DYNAMIC PROPERTIES TEST 


The prlmai’y purpose of tlie propellant test was to determine real moduli 
over a range of dynamic parameters anticipated for the MVGVT and the Orblter 
Flight Test (OFT) vehicle service conditions. Pertinent test considerations 
wei’e noted and resulting ' lata were combined and smoothed to conform to 
specific trends. 


A. Test Considerations 

Because of the urgency for supplemental propellant data, general teat 
requirements are baaed on relevant published Information presented In Refer- 
ences 1 through 5 and specific test consideration for facilities available follow. 

1. Specimen Design. The lap shear type used by Rockwell Inter- 
national [2] appeared to best distribute applied test loads uniformly. It had the 
added advantage that Its tester could accept a selection of large electrodynamic 
shakers to cover the desired range of dynamic strain amplitudes. 
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The specimen size had to be large enough to represent tlie properties 
of the bulk grain and had to repress effects of nontyplcal variations that might 
be manufactured Into the grain by flow and settllnp; patterns during casting and 
by localized temperature gradients during cure. Ihe shear area could not 
exceed a 2- by 4-ln. size because of the limited shaker force available for 
testing at the required low temperature and strain amplitude. The 1-ln. 
thickness used by A. J. Yorgladls [1] seem.ed to satisfy the Instrumentation 
requirements and minimized the bonding material property affects over the 
propellant properties. 

2. Static Strain Simulation. Usually tost specimen are cut out of 
batch samples cured In flexible pasteboard cartons which preclude cure and 
thermal shrinkage strains. Consequently, their properties ai’a somewhat 
different from propellants cured and stored In rigid motor casings. 

As the cast propellant cures and drops to ambient temperature, T, the 
grain shrinks linearly bj 

[1 +0.055(T - 146)] X 10“3 (16) 

toward the bonded casing, and the Inner radius Increases and creates radial 
and circumferential strains 


= H(1 + v) j^l + 

and 

respectively. The Poisson’ s ratio is taken as 0. 5. Combining these strains 
and averaging over the grain cross section having an outer radius twice the 
Inner radius results in a static strain of 


y c- 6H 
's 


(17) 
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It was recognized that this static strain develops over many hours which 
allows the stress and, tliorofore, the associated modulus to decay nearly as 
rapidly as It is applied. Also, the static strain then remains relatively con- 
stant through the SRB storage life with a static modulus less than half the 
Initial value. The significance of the static strain Is that, once Incurred, the 
material stinicture may be permanently weakened to the extent that subsequently 
imposed strains will exlilblt a lower modulus due to separation of filler-binder 
adhesion of the uniform mlxtui’o. 

To simulate static strain on test specimens at a real time rate Is 
unthinkable for the production type testing required. Instead, a parametric test 
was requested In which several levels of static strains were specified to be man- 
ually Imposed and the modulus allowed to relax a minimum of 20 minutes prior 
to dynamic testing. If the modulus did not appreciably change at static strains 
up to 3 percent, this parameter would be Ignored, 


B. Test Data 

Two distinct sources of data wore used to compose the required dynamic 
properties of SRB propellant. Thlokol provided the statistical mean properties 
using one strain amplitude while R1 provided supplemental data at lower fre- 
quencies and over a range of dynamic shear strains as shown In Table 1. 


TABLE 1. PROPELLANT DATA SOURCE 



Thlokol 

Rockwell 

International 

Reference 

2 and 6 

7 

Type Specimen 

Gottenberg 

Lap Shear 

Number of Specimens 
per Temperature 

18 or more 

1 or 2 

Frequency Range 

5 to 50 Hz 

1 to 20 Hz 

Dynamic Strain Amplitude 

0.008 

0.001 to 0.014 
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1. Combining Data . Combining the two sources of data consisted of 
using the real modulus mean value at 5 Ha and 0. 008 dynamic strain provided 
by Thlokol and superposing the 0. 001 to 0.014 spread of dynamic strain pro- 
vided by RI with both coinciding at the 0, 008 strain curve as presented In 
Figure 4. 

The real modulus obtained at 5 Ha from RI data was within the deviation 
range of Thlokol data except that the RI Inert propellant was slightly higher. 

The real modulus deviation from specimen to specimen and batch to batch was 
approximately 10 to 20 percent with tlio higher temperature propellant demon- 
strating the larger deviation. In general, tlio Inert and live properties were 
very similar. 

The dynamic strain slopes of both sources of data were almost Identical 
at 5 Hz. Spacing between dynamic Btral!’i Increments was noted to be logarithmic 
as expected of rubberllke materials. This was also noted In Reference 1. 
Dynamic strains were plotted smooth to compensate for spurious data points 
of the limited number of RI specimen used, but patterns and specific trends 
observed were preserved. 

2. Static Strain . The effect of static strains to simulate cure and 
thermal shrinkage was not Incorporated In Flgui’e 4 and Is not recommended for 
future analysis imtll It can be supported by a detailed test program. The real 
modulus was expected to decrease with Increasing static strain [5] but the 
reverse occurred. This could be Indicative of unaccounted compression which 
Increased the friction between the binder and filler, probably because 

of the specimen design and loading method. The preload application and relax- 
ation rate needs to be developed for production typo testing to reasonably simu- 
late shrinkage conditions. 

3. Data Results . The data provided satisfied all requirements except 
for static strain, which Is a refinement generally Ignored for probably the same 
difficulties encountered here. The composite development was bettor tlian 
expected because the two sources of data were reasonably compatible at the 
common superposition point of 5 Hz and 0.008 dynamic strain. Though dynamic 
strain properties only for 76“ F ambient are presented for Inert propellant, there 
are enough similarities between Inert and live properties to admit sufficient 
extrapolation and trends to estimate properties at any temperature within the 
40“ and 90 “F range. 
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DYNAMIC SHEAR MODULUS. G X 10 ^ PSI 











V. GRAIN MODULUS DETERMINATION 


Having Identified and Interrelated tlie parameters Influencing the gi’aln 
dynamic real modulus, their utilization and procedure for determining the 
modulus Is now demonstrated In application to MVGVT and OFT vehicles. The 
analysis Is limited to the aft attach ring behavior and Is consistent with the 
study scope. The MVGVT Is treated first because of the required teat data to 
define tlie actual dynamic behavior of the Shuttle mated system and resulting 
environment on the SRB at the ring station 1517, The lift-off response was 
then assumed to be dominated by the MVGVT system fundamental mode. 


A. MVGVT Vehicle 

Li tracking those parameters directly affected by teat environments, 
two MVGVT configurations were noted; a 0-pal and a 600-psl combustion 
pressure simulation. Because the two configurations were tested a month 
apart, temperature parameters were expected to bo significantly different. 

It was also noted that two fundamental modes were excited. The sym- 
metrical mode consisted of the SRB elements vibrating with opposite angular 
accelerations at the aft attach ring level. The nonsymmetrlcal mode had both 
SRB vibrating with similar angular acceleration. Therefore, the pitch exci- 
tation parameter, g^^, of the two modes were different. 

Figure 5 summarlzeis relevant MVGVT test data from which environment 
parameters of grain temperature and pitch excitation may bo determined. 

Configuration 1, with no combustion pressure simulation, was tested 
on October 26 and 31. For three consecutive days prior to the test, the average 
temperature was 65® F which was also believed to be that of the propellant grain. 
The second configuration, wltli a 600-psl pressure simulation, was tested on 
December 1 and 2 with a probable gx’aln temperature of 50* F which was also 
based on a 3-day average ambient; Propellant properties at 50® and 65® F were 
constructed and are presented In Figure 6. 

The average pitch excitation parameter, acting on the SRB ring 

along the Z~axls Is bound by the acceleration diagram shown in Figure 5 
and Is calculated by 
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Figure 6. H-18 propellant properties extrapolated to MVGVT grain temperature. 
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" II. 



whoro 


/j - 0. r>(278Z + -102Z) ^ 

(jjHj - 0. 5(278Z - »02Z) ) 

- 0. r>(277Y + 278Y) 


( 1 «) 


(ly) 


and thu alphatuimerics In paronthuaes uru acceleration data of dcalKTiatc^ 
ai-celei'ometei 3 ( l•■l^ture fj) . 

UesultH of parameters from MVOVT lest data are aummari/.oci In 
I al)le 2 where is the maitni system fundamental frequency. This frequency 

Is plotted as a fuiu tion of jillch excitation, in Fljture 7. Sul)scripts M 

refer to matcxi vehicle conditions. 


TAHl K 2. MVC.VT KN VHtON MKN I A I l>Alb\MKTKILS 


(’otifljfiiration 

Simulatixl 

Pressure 

(psl) 

tlrain 

Temperature 

CF) 

— 

Test 

M(xie 

Symmetrical 

System 

(j) 

M 

(Hz) 

SUH 

Excitation 

1 

0 

dfj 

■ 

Yes 

2.059 

0.090 




« 

No 

2.29.5 

0.011 

') 

<)00 

f)0 

•17 

Yes 

2.o:$9 

0.099 





-18 

No 

2.215 

0.090 
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Figure 7. MVGVT propellant grain modulus. 


The next step was to determine the frequency and modulus^ G' » of a 
free vibrating ring-grain structure if excited by the acceleration experienced 
in the MVGVT of Table 2. Using the modal analysis program*, two arbitrary 
moduli were selected to establish curves (a) and (b) of Figure 7 which extends 
through the excitations of MVGVT. Table 3 presents output data required to 

calculate P of Equation (13) as required by Equation (11). Table 4 presents 
the output data from the static analysis used to determine the metallic ring 
stiffness, C, of Equation (14) and the static constant, D, of Equation /"12). 

Both are also required by Equation (11) . The dynamic strain \vas selected 
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TABLE 3. MODAL OUl'PUT 


Simulated 

Presaure 


G 

Assumed 

b> 

(Ha) 

T3 

(209) 

n 

F4 

(Eq. 8) 

P 1 

(Eq. 13) 

415 

1.40 

0.960 

0.01149 

1.915 X10"‘ 

800 

1.58 

0. 962 

0.01139 

1.172 X 10* •' 

415 

1.87 

0.901 

0.01033 

2.951 X 10"* 

800 

1.92 

0.878 

0.01008 

2.688 X 10“* 


TABLE 4. STATIC SYSTEM DATA 


G 

Asflumed 


G -0 
s 


G^ » 400 
0 


G » 0 
s 


Output 

Deflection 

(in.) 

C 

(Eq. 14) 

1 

D 

(Eq. 12) 

AZ -16.58 
s 

210 

4920 

AZ^ - 5. 71 



AZ -2.82 
8 

AZq = 2.29 

3470 

13,830 


2 ^ 


























from the MVGVT propellant data of Figure 6 for each set of modulus and natural 
frequency presented In Table 3 at appropriate temperatures obtained from 
Table Jj. These proportlea are summarized in Table 5, which includes tlie exci- 
tation calculated by Equation (11) to induce the associated dynamic strains. 


TABLE 5. MVGVT FREE VIBRATION SYSTEM 


Configuration 

Temporaturo 

(^r) 

G' 

0) 

^d 

(Eq. G) 

g 

(Eq. 11) 

Curve 

1 

65 

415 

1.40 

0.025 

0.265 

a 



800 

1.58 

0.002 

0.035 


2 

50 

415 

1.87 

0.20 

4.90 

b 



800 

1.92 

0.02 

0.538 



The data presented In Table 5 are all tliat are necessary to establish 
curves (a) and (b) for moduli and frequencies over a range of excitations, g, 
plotted In Figure 7. It Is also of Interest to note that all parameters of any 
free vibrating system with a known Initial excitation g cannot ho derived by work' 
Ing backward through the established equations presented In paragraph m.B. (2) 
and the propellant properties of Figure 6 without the aid of curve (a) or (b) . 

The ring-grain structure In the MVGVT vibrates with the vehicle system 
fundamental frequency, which may not be the same as the free vibrating system 
though selected wltli the same excitation as In Figure 7. Therefore, the moduli 
are not the same. To determine the grain modulus acting during MVGVT, the 
relationship of motion between L.e two systems Is investigated. The subscript 
M will refer to tlie forced system in the mated vehicle as before and n will 
refer to the free vibrating system. 

The displacement of the free system witli Initial conditions set by g^ Is 
calculated from Equation (9) as 
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2 = tjT3 ^ (C +G» ) T - : — 

n ' ' 0' (C + G|^) 


or 


D 


^n“ ®M n/Y (C + G») 


The displacement and acceleration relationship of the forced harmonic 
system Is obtained In the usual way of differentiating displacement twice so that 


_ A._ (32.2 X 12) 


Z. 


M o>' 


M 


(27t)‘ 


Solving for g^^ In the displacement equation of both systems and equating, the 
ratio of the two systems Is obtained as 


_ (^/2 X 32.2 X 12) 

^ Z Do)^ (27r)^ 

n M 


( 20 ) 


From curves (a) and (b), the free vibration system modulus and fre- 
quency Is obtained for the g^. Using these two properties In Figure 6, the 

dynamic strain Is obtained. Since the dynamic strain of the propellant grain Is 
proportional to displacement, the dynamic strain of the forced system Is 
calculated by 


Using this dynamic strain, y^, and the frequency of the MVGVT system, 

ct>M» Figure 6, the desired modulus G^^ Is resolved. Data of the two systems 

used to resolve the MVGVT modulus for each mode are presented In Table 6 and 
moduli, are plotted In Figure 7. 
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TABLE i). MVGVT PROPELLANT GRAIN DATA 



Mode 

"n 

G’ 

n 

7 

'^n 


r 

'^M 

G' 

M 

0,030 

6 

1.6 

860 

0,0017 

2,039 

0.70 

0,0012 

1000 

0.011 

8 

1.7 

1200 

0.0004 

2.235 

0.79 

0.0003 

1400 

0,039 

47 

2.0 

1700 

0.0016 

2.039 

1.24 

0. 0018 

1600 

0,030 

48 

2.0 

1800 

0.001 

2.235 

1.07 

0.0011 

1900 


B. OFT Vehicles 

Unlike tlio MVGVT, the OFT flight environments at lift-off are defined 
by b I'D ad ranges, but combustion pressure of the SHB may be assumed to be 
600 psl at maximum excitation so that only one configuration will be analyzed. 
The extreme ambient temperatures specified are 40* and 90“ F. Two funda- 
mental modes for the mated vehicle system that Influenced the pitch excitation 
parameter, g^^, were noted to exist: symmetrical and nonsymmetrical. Both 
modes will be considered for the OFT flight environment. Because ground 
winds and engine thrust at lift-off are variable contributors to pitch excitation 
of the mated vehicle, a range over one oixlor of magnitude of acceleration will 
be Investigated. Consequently, the OFT g^aln modulus In tlie aft attach ring 
will be resolved parametrically. 

The pitch excitation pai'ameter of Equation (18) may be more meaning- 
fully defined for flight vehicles by referring it to the anticipated Z -acceleration 
at the attach fitting which corresponds to the accelerometer number 278Z 
location on the MVGVT. It Is an Interface point common to the Shuttle elements 
and design loads community. Then, OFT pitch acceleration acting on the SRB 
center is 

= , ( 22 ) 


where g is the anticipated acceleration at the fitting. Substituting MVGVT 
values for tlie symmetrical and nonsymmetrical modes, one obtains 
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0.95 


K 


B 


M 0.039 




1 0.041 

*278Z 


g 


B. - 


M 0.030 


“ *278Z 


0.81 


( 23 ) 


respectively. 

Table 7 summarizes the lift-off environment and parameters to be 
applied to the OFT modulus determination that covers extreme yearly opera- 
tional conditions. 


TABLE 7. OFT ENVIRONMENT PARAMETERS 


Mode 

Grain 

Temperature 

(•F) 

System 

U) 

M 

B 

(Eq. 22) 

^F 

Symmetrical 

40 

2.039 

0.95 

0.1 to 1.0 


90 

2.039 

0.95 

0. 1 to 1.0 

Nonsymmetrlcal 

40 

2.215 

0.81 

0. 1 to 1.0 


90 

2.215 

0.81 

0. 1 to 1. 0 


To establish curves (a) and (b) in Figure 8 of a free \1bratlng ring- 

grain structure as was done on MVGVT, Tables 3 and 4 were applicable again. 

Figure 4 was used to select the live propellant properties of y . listed in 

d 

Table 8. The excitation, g, was calculated from Equation (11). Curves (a) 
and (b) of moduli and frequencies were extended through the excitation, g, 
range of 0. 1 to 1. 0 in Figure 8. 
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TABLE 8. OFT FREE VffiRATION SYSTEM 


Temperature 

(*F) 

G* 

fa) 

^d 

(Fig. 4A) 

m 

Curve 

40 

415 

1.87 

0.28 

6.86 

b 


800 

1.92 

0.028 

0.76 


90 

415 

1.87 

0.0012 

0.03 

a 


800 

1.92 

0.00012 



0.003 



Finally, the excitation acting on the ring center was calculated by 

Equation (22) for select anticipated excitations at the fitting. Free vibration 

moduli and frequencies were inteipreted from Figure 8. The modifled dynamic 

strain for the OFT grain was calculated by Equation (21). Propellant grain 

data of the two systems are presented in Table 9. The desired OFT moduli, 

G’ , at 40* and 90* F are plotted in Figure 8, which is almost equal to that of 
M 

the free vibrating system because the frequencies are almost equal. 


VI. CONCLUSIONS 

The objective to determine the grain dynamic modulus in the SRB aft 
attach ring during lift-off was accomplished through the development of a com- 
prehensive analysis which considered the variation of dynamic strain property 
of propellant. Results of a study for fli^t vehicles are plotted in Figure 8, and 
pertinent observations were as follows: 

1) The most significant variation of modulus is due to temperature. 

The real modulus increases by a factor of five from summer to winter launch. 

2) Increasing the pitch excitation at lift-off due to increasing ground 
winds and engine thrust decreases the real modulus by less than 60 percent in 
the range of Interest. 
















TABLE 9. Ot'T PHOPELLANT CHAIN DATA 



B 

Si 

Temperature 

(•F) 

u 

n 

B 

^n 


B 

Si 

C 

M 


0,95 

0.095 

40 

2.0 


0.0032 

2.039 

^9 

0.0038 

1450 

1 



90 


300 

0.0045 


0*91 

0.0040 

300 

1 

0.81 

0.081 

40 


1650 

0.0029 

2.215 

1.04 

0.0030 

1650 

1 



90 


310 

0.0048 


0.77 

0.0037 

330 


0.95 

0.285 

40 

1.95 

1100 

0.01 

2.039 

1.10 

0.011 

1050 




90 


210 

0.02 


0.88 

0.018 

210 


0.81 

0.243 

40 


1150 

0.01 

2.215 

0.94 

0.0094 

1150 




90 

1 

t 

220 

0.019 


0.75 

0.014 

240 

1.0 

0.95 

0.95 

40 

1.90 

1 

7.50 

0.034 

2.039 

1.01 

0.0345 

750 




90 

1 

1.50 

0.075 


0.87 

0.065 

160 


0.81 

0.81 

40 

! 

1 

800 

0.032 

2.215 

0.87 

0.028 

850 




90 

1 

t 

1 

* ■■ - 

160 

0.070 


0.74 

0.052 

170 


3) The highest pitch excitations at the fitting anticipated for any launch 

day is 0.8 g and the lowest with negligible ground winds is 0.4 g . Design 
M M 

modulus for any day launch should range from 170 to 1000 psi. 

4) SRD combustion pressure of 000 psl was used which Indirectly doubled 
the grain modulus. Pressure buildup of less than 600 psl during lift-off will 
decrease the structure stiffness. 
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5) Grain dynamic strains resulting In the aft attach ring range from 

0. 01 to 0. 07 for flight vehicles. Dynamic strains forward of the ring are much 
less. Dynamic properties data from propellant contractors should Include a 
broad range of dynamic strain tests rather than the present practice of a 
single amplitude provided by References 2 and 6. 

6) Simulation of SRB cure and stoi’age static stralr.s on teat specimen 
should bo formulated for production type testing. 
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APPENDIX 


DYNAMIC STRAIN COMPUTER PROGRAM 
FOR SRB GRAIN CROSS SECTION 

Tho strain analysis was programmed In Basic Language for the 4051 
TEKTRONIX computer. Symbols used are the same as defined In paragraph 
m.B. (1) of this report %vlth tho exception of: (1) Y and Z replace the T2 and 
T3 elgenvectox’s used In Equation (l)} (2) K la an optional conversion factor 
wltli no present assignment; and (3) P(I) andP(M) replace angles ^{l) and 
^(2) of Equation (l) and Figure (3). 


ZC 



FIND 1 

OLD 

LIST 

18® RFItf'.SK Propellant Grain Dvitanlc Strain Proprtxm 
161 DCLETE Y,2,A>B»R,8 

192 SET DEGREES 

193 UIEUPORT 29,119,9,99 

194 Din Y<3,17>,2<3,17>|A(2>,B(3>,R<3>,8<2> 

195 DATA 1,8.495,16,2,39.917,48,72.74 

196 READ K.U.H.HfR 
187 FOR J*1 TC M 


les 

18!» 

118 

111 

112 

113 

114 
* 115 

1 16 
! I"* 
ltd 
t 


126 

121 

122 

123 

124 

125 
_ 126 

•^123 

129 

T,4i4M lie 

131 


S< J>*R< J»l)-R<J) 

A* J^=PIt<R<J♦l>t2-R(J)t2)/N 
B< J>*2tPUR< J)/N 
NEXT J 

6vHfl)=?tPI4R<n4l)/M 
PRINT *L 
PRINT "It- 

PRINT "Z Input nunber of sectors 
PRINT "f * iiuMber of annular rinns 


Profran Hemi 1. Input Run "l 

. *1 

Defouit * 16_ 3. Inn"; 
DefoclT » 2. 4. Ir.pu* disp*l 
fault * 1_ 5. P';iSiC'''r r*. 


FRIMT O’ j Dofi.^lt » .495^ 6. Irou* r;nj ro*; 

^RIliI Default - 30.917, 48, 7?. 74^ 7. f-put dispia’; 

PRINT •'v'tKjTits fron tape file _ 8. Input displacenents”! 

PRINT " ’on keyboard. 9. Output displactficntt to tape file.18. Ou*l 
PRINT "tr'it elenent nofieric1clure.ll. Output calculated clencnt fco*! 
PRINT *n«'try.l2. Output propellant dynanic strains .13. Exit pree*! 

PRINT ”ran Inpat line naiiDer of desired proyran resp on se : *| 

INPUT R1 

GO TO R1 OF l»,132»137,169»168,194f229,212,238,142,194,244»127 


PRINT '..Run ID • "I 
INPUT TS 
CO TO 113 



(-0 






R's 




32 

33 

34 

35 

36 

37 

38 

39 
48 

41 

42 

43 

44 

45 

46 

47 

48 

49 
58 
51 
5? 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 


PRINT • Hmil»cr of sectors ■ *1 

INPUT N 
DELETE Y,2 

DIN Y(N^1,N41>,2CN^1,N^1> 

CO TO 187 

PRINT ■ Hunber of annular rints ■ 

INPUT N 

DELETE Y Z A B R S 

DIN Y<H+l,N^l)Iz<N4l,N^l),A<f1>,BfN^l >.R<H^1),S<H) 
GO TO 154 

PRINT "LRun ID * 'iTl;" Element Nomenclature" 

MINOOH -R<n4l),R<H4l),-RfH4l),R<«4l > 

FOP 1=0 TO N-1 
COSUB 171 
NEXT I 

PRINT 832,21:115,45 

PRINT 832,20:112,47,112.43,115,45 

PRINT ■ Y" 

PRINT 832,21:65,95 

PRINT 832,20:63.92,67,92.65,95 

PRINT "tZ" 

GO TO 209 
P jP 1=1 TO «♦! 

IP I /I THEN 158 

PRINT • Inner radius* 

CO TO 162 

IF I*N^1 THEN 161 

PRINT "Intermediate radius* "| 

CO TO 162 

PRINT "Outer rodius* "I 
INPUT R(I) 

NEXT I 
CO TO 187 

PRINT " Oieelecement Scale, K* "| 

INPUT K 




_1€7 GO TO 113 

1€8 PFIHT • Input Poisson's ratioi *1 

169 INPUT U 
_170 CO TO 113 

iv, 171 REKARK Sector Plotting Subroutine 
VrmM J72 Pl«36esi^H 
^*•^171 P2-Pl436e/^N 

174 NOME R<l)tCOS<Pl),R<l)SSIN(Pl) 

175 DRAM R<«4l)tC0S<Pl),R<H4l>SSIN<Pl> 

176 NOUE R<l)tC0S<P2)|R(l)tSIN(P2) 

177 DRAM R<H^l)tC0S<P2),R<H4l)*SlN(P2) 

178 FOR J=1 TO H 

179 nOUE R(J)*COS<Pl>,R(J>tSINfPl) 

180 ORAU RU>tC0SfP2),R(J)tSIM<P2) 

181 R2=<R(J)^R(J+l))tCOSU80/H>^2 

182 NOUE R2tC0S<<Pl+P2)/2>,R2»SIN<<Pl4P2>/2) 

183 PRINT -E-U^Htl 

184 NEXT J 

185 NOUE R<N^l>tCOS(Pl>,P(H^l)tSIN^Pn 

186 PRAM R<N+ritC0S(P2>,PfM*l >tSIN(P2 • 

187 RETURN 

REHARK Define Sector Uindous 
WirMiM 189 Rl*(R^H4l)4RM>)#COS«'360*vI^0.5^/M - '2 

190 R2*(RvH4l)4R<l))tSIN(360t(I^0.3>N»'2 

191 F*R<H+1)-Rf 1) 

192 UINDOU R1-F,RUF,R2-F,R2^F 

193 RETURN 

(U*p«v^94 PRINT *LPun ID • ■|T$|" Calculetcd Elcnent Geometry* 

(m) 195 1*1 

196 G0SU6 188 

197 GOSUe 171 

198 FOR J-1 TO H 

199 NOUE R(J>tC0S(S48^N>tC08(188/N>,R<J>t8INCS4SrH>tC0S(188/N> 
2M PRINT -B'lJf- -IBCJ) 

2fl R2-<R(J41>^R(J>>/2 






292 

294 

29 Z 

29 € 

29 ? 

296 

299 

219 

.211 

2\2 

213 

214 

215 

216 

217 

218 
.219 
2?9 
221 
222 
223 


224 

225 

226 
227 

_228 

229 


3a^c 230 

03 231 
^ 232 

233 

234 
^35 

236 


nOUC R2tCOS<3<0/N>,R2tSIN(360^N> 

PRINT "S'lJl** -ISCJ) 

NOUE R2tCOS(360/N>,R2tSIH(360/^N>-0.02)t<R<N4>l>^R<l>> 

PRINT •A*|J|*» -lACJ) 

NEXT J 

NOME P(H^l)tCOS(54e^N)tCOS(18e/N>,R(N^l>tSlN<S4e/N>tCOS(180/N> 
PRINT •B"lN^l|-« •|B<N^1) 

PRINT »32,2l:iie,5 
INPUT P$ 

GO TO 113 
FOR 1-0 TO N-1 

PRINT 'LRun ID * "iTIi" OisplaccMcnt Data” 

GOSUB 188 
GOSUP 171 
FOR J=1 TO H+1 

NOME R< J>tCOS(Pl),R(J)»SIN<Pl) 

PRINT ’dY- •; 

INPUT Y<J,l4l) 

NOUE P- J‘*tCOS(Pl>,R( J)tSINfPn-i5..>2t(R<H4l WRf n) 

PRINT ”dZ» 

INPUT ZfJ.I^D 
IF I>1 THEN 226 
Y<J,N^1 )«Y< J, 1 ) 

2(J,N^1)-Z<J, 1) 

NEXT J 
NEXT I 
CO TO 113 

PRINT ■ Input file nunber - ”| 

INPUT F 

IF F<«0 TNEN 229 
FIND F41 
INPUT 033: Y»Z 

PRINT ”LdY'* iron File i*|F,Y 
INPUT P$ 

PRINT ”LdZ'f iron Flit •"IF.Z 


*r 3 

is 




X 


aTc 237 CO TO 2tS 

230 PRINT * Output flit nunbcr ■ "I 

7^ 239 INPUT F 

^ ' 248 IF F<«C THEN 238 

_241 FIND f*l 
242 PRINT 833: Y|Z 
? -3 CO TO 113 
i.^.4 P-8 

245 F2-8 

246 F4«8 

247 S2»8 

248 PRTMT "LRun ID • "lT$i'__Ele»e«t Ul U2 *1 

249 PRINT • U7 H F3_* 

258 FOR I-P TO N-1 

251 FOR J-1 TO H 

252 PRINT •E"lJ4HtU'^K* 

253 Pl-368tIxN 

254 P2-PU368/N 

255 Ll«V<j4l,I^l)tC0S<Pl)4Z<J*M,l4l)tSIN<Pl) 

256 L2*Y<J4l,I>2)tC0S<P2)^Z<J4l, I^2)tSIN(P2) 

257 L3=Y<J, I^l)tCOS<ri)+Z<J,I+l)»SIN<J»l) 

258 L4*Y< J, l42)tC0S(P2>^Z< J, l42>tSIN<P2) 

259 Hl»Y<j4l,l4l)tSIN<Pl)-Z<J4l,I+l)8C0S(Pl) 

268 H2-Y<j4l,I^2)«SIN<P2)-ZOl,I>2>tC0S^P2) 

261 H3*Y<J,l4l)tSIM<Pl)-Z<J,l4l)tC0S<Pl) 

2'2 «4-Y< J. l42)tSIN<P2> -Z< J, l42>tC0S(P2) 

263 U1*<L14L2-L3-L4»'2/8<J) 

264 U2«<(mHI2>/B<J4i)40l3HM>/9(J>)/2 

265 U3«((lll4ll2-ll3*N4>/8(J>4ai-L2>^B(J4l>4(L3-L4>/e(J>>/2 

266 Cl-(Ul4U2>/2 

267 C2-8QR<(U1-U2>t244l3t2»^ 

268 Fl-Cl4C2t8GN<Cl> 

269 F2*F248M<F1>8A<J> 

278 F3»Ci-C2tSCN(Cl> . 

271 F4-F44(Flt24F3t2>SA(J> 



272 

273 

274 

275 
27€ 

277 

278 

279 
288 
281 
282 

283 

284 
2r,5 
286 

287 

288 
289 
298 

291 

292 

293 


S2*S24fitA<J> 

PtlNT USING 274:Ul,U2,U3iFl|F3 
INAGE 11D.8D»4<4D.80> 

IF J^NtI«24 THEN 289 
NEXT J 
NEXT I 

F2«F2^I^(R<IH1 )t2-R< 1 )t2) 
F4*KtSQR(F4/PI/’<R<«4l )t2-R( 1 )t2) ) 
S2*S2/PI^<R<H4l)t2-R<l)t2> 

P$**N€t tension* 

IF S2>8 THEN 285 
Pf«*H«t conpression” 

IH8GE L,*F2, F4 • *»20.8D,4D.80, 14X, 15A/47X, 

PRINT USING 286:K,V 

IHPGE L,*X, U » *,4D.8D,4D.8D 

PRINT USING 284:F2,F4,Pf,S2 

CO TO 289 


P«I4N-1 

PRINT USING 286:K,U 
INPUT P$ 

GO TO 248 
END 


■<*,20.80, ->- 


'Hi" '•^l-' 
: \-d T\ K1! ' 


O 


Pro9rQn Mer^u 


1. Input Run ID 

2. Input number of sectors Default ■ 16 

3. Input number of annular rings Default * 2 

4. Input displacement scale Default ■ 1 

5. Input Poisson's ratio Default • .495 

6. Input annular ring radii Default ■ 3e.917f 48i 72.74 

7 . Input displacements from tape file 

8. Input displacements from Kewboard 

9. Output displacements to tape file 

10. Output element nomenclature 

11. Output calculated element geometry 

12. Output propellant dynamic strains 
13* Exit program 

Input line number of desired program response 1 1 


Run ID « FILE 
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